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CALCULATING  AN  ACOUSTIC  FIELD  IN  LAMINAR-HETEROGENEOUS  MEDIA 
V.  M.  Kudryashov 

The  practice  of  calculating  an  acoustic  field  shows  that  the 
better  the  components  (whose  sum  describes  the  field  to  be 
calculated)  correspond  to  the  spatial  structure  of  the  field,  the 
more  effective  the  selected  means  of  calculation. 

Let  us  examine  a  field-calculation  method  convenient  in  the 
region  where  the  zonality  in  the  field  structure  applicable  to 
uniaxial  channels  appeared.  For  this  let  us  introduce  a  cylindrical 
coordinate  system  (r,  0,  z)  (Fig.  1)  and  examine  the  acoustic  field 
in  a  layer  big  £  z  <  H  filled  with  a  medium  which  has  no  shear 
elasticity,  in  which  the  square  of  the  refraction  index  n2(z)  and 
its  derivatives  of  an  order  no  lower  than  the  first  are  continuous 
functions  of  the  z  coordinate.  In  the  half-space  z  >  H  we  set 
n  (z)  =  n  (H)  +  2a(z  -  H),  n(H  -  0)  5  n(H  +  0).  The  acoustic 
potential  (of  velocities)  Y(r,  z,  zQ)  exp(-i2nfQt)  should  satisfy 
the  wave  equation  everywhere  except  point  r  =  0,  z  »  zQ  in  which 
there  Is  located  the  point  source  of  a  sound  operating  under 
stationary  conditions  at  frequency  fQ  >>  C(z)  |  (d/dz,!h2(z)  |  ;  In  the 
vicinity  of  the  source 

V(r.  in)  exp(lk„L)IL, 

where 

jTp  _  o. 
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In  addition,  lim  Y(r,  z,  zQ)  ■*  0.  The  medium  is  bounded  from 
above  by  the  plane  z  ■  hg  on  which  ¥(r,  hQ,  zQ)  «  0. 

In  this  case,  according  to  L.  M.  Brekhovskikh  [1], 

/«*2  “t‘* 

V(r,  t,  (<!',(*„;  «)  -/>(«) <!>,(*„;  *)|  X 

X I'M*  *)-f  («)*,(*:  «)||i  - /5(«)?(»)p' /yj" (*/»)«</*, 

where  ^(z;  k)  and  4>2(z;  k )  are  two  linearly  independent  solutions 
of  the  equation 

jp  4»  (*)  +  kl  [«*(«)  -  »'] (*)  -  0.  ( 2 ) 

Magnitudes  q(<)  and  p(tc)  are  determined  from  the  boundary 
conditions  on  the  upper  and  lower  boundaries  of  the  medium, 
respectively,  or  from  radiation  conditions  if  the  medium  is  unbounded. 
In  the  case  in  question  we  get 

?<«>  «)/•!»,  (I>„;  «): 


•  «,  («;  M  -  !L  [«|!{  (-i„)///J»  (,„)]  |nJ  (II) ...  ,J|-'  » 

1) 

dF*i  (*;  *) 

L-w 

*»<//;  »>  7- l»l!i (!»)/" •,i(iw)|v,i»J(//) 

IF  *1  (*>  *> 

-*»r- 


The  integration  contour  in  equation  (1)  circumvents  those 
positions  from  below  in  which  Re  k  >  0,  Im  k  >  0,  and  circumvents  the 
remaining  singular  points  of  a  subintegral  function  from  above. 

If  z  <  h,  (d/dz)n2(z)  >  0,  then 

M*  «WW  Kil«*W—Twe*p[<—  <-&-)<■  +  A,). 

(3) 

*•  (*:  *> "  /  J  W  “  *’•  (*»  -  <tt)0  +  A,)  - 

-/"<'>  e*P  (1  +  A,)] , 
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1>, (*;  «)-/  ■'  («)  yifc  |«»(*>  -**)-u,exp(-  <»  +  <  -n-)(1  +  A,). 

-  -Tp<afg»<-y-: 

(3) 

«*,(*;  »)=  |«JW  -  *!r 1,4 [/(«)-'  exp  (-<»»  +  <  TJ7  0  +  A«)  ~ 

—  /(«)exp  (Ao  —  +  A,)j ,  -j-<  arg«<-j-  ■ 

In  formulas  (3)  the  designations 

t 

»  =  »(*;  ,)  =  *„  J  /«’(«) 

#1 

are  used,  where  is  the  root  of  the  equation  n(z)  »  k,  which 
satisfies  the  requirement  z^  <_  h,  if 

Im  •=  0;  ilmZ|  =  llm  (*  —  •)■ 

*•-*  *•  (A)  I  -*0 

Since  oi ( z ;  <)  is  a  multivalued  function  of  k  with  branching  points 
k  =  ±h(z)  ,  we  assume 

nrgu^O,  at  «a  <  «*(*)  <  1,  lm«J «  0; 

at  «(*)<*<!,  lm»=.||in(0-»); 

•-♦0 

nrgm=  —  -it,  at  «(*)<*<  1,  lm«  —  |lm(0+»). 

•-•0 


The  whole  picture  is  symmetrical  relative  to  the  point  <  =  0 . 

The  introduction  of  factor  f(<)  into  formulas  (3)  is  associated 

2 

with  the  presence  of  a  maximum  n  (z)  at  z  =  h;  f ( < )  assures  a 
regular  behavior  of  the  right  sides  of  expressions  (3)  on  the  entire 
plane  k. 


where  oi0  =  to(z2;  k);  z0  Is  the  second  root  of  the  equation  n(z)  * 

=  k  z2  >  h  at  Im  z2  =  0;  f(x)  ^  1  if  |w2|  »  1,  |  arg  u>2 1  <  it.  And, 
finally , 

A, («)-»(--);  |A,(*)K<  l.  |»|>l. 
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Asymptotic  expressions  for  p(k)  have  the  form: 


/»(«)=  V  (.)/-*  (*)exp  [/2K  -  «„)  -  l  -  *-] . 

( ^a) 

p  («)  -  V  (*)  exp  (/2U.J  -  l t)  [‘  "  /  2  <*)  e*P  ( 
*  .  3* 

-<2«„+^)j, 

(4b) 

3n 

-5- , 


where 


(•+*)' 

,  f,  1  l  /*w«p («»„  —  / t) 

p  1— /*<»)««p(/2»w  — I-J-)  _ 

1 +<-)['- 

/’MexP(/2«„-/-f)]j  '[>+°(-^ 

p  is  the  density  of  the  medium  at  hQ  <_  z  <  H;  pg  is  the  density  of 
the  medium  at 


*>//;  UH  =-«(«;«): 


“(*:  «)  =  *oj  V ft’ ({)-»*(«; 


arg  u(z,  H)  on  the  upper  sheet  of  the  Riemann  surface  is  set  by  the 
same  conditions  as  for  <d(z;  kt).  In  the  derivation  of  formulas  ('4a) 
and  (<4b)  it  was  assumed  that  2b  »  |g-zn2(z)|. 

It  is  easy  to  see  that  the  right  side  of  equation  (I)  is 
identically  equal  to  the  right  side  of  the  equality 

w.=i 

V-=<r,  *,  *,)-  2  *<  *>)  +  ^’>Ar  *.  *«.).  (5) 

where 

4  •  +  !• 

t,  J  F„{r,  z,  ( 6a) 


V*.(r,  *.  «„) 


+Vh 

J  F*,(r,  z,  z„-  .)  1 1  -  p («) q Ml- rf«. 

--  +  !« 
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Here 


F„{r,  z,  •)  =■  -j- •».(«•;  «)|®,(«; «)  «)|  //•(*/"): 

F„(r,  t,  «)  =  _/  |<J>,(*0;  «)  -  q (*) «J>,  (*,;  «))  X  (  7 ) 

X  («>,(*; «) -*(«)*,<*:  «)f*‘ V'MHff’flW.  N>  1. 

The  representation  of  expression  (5)  is  very  convenient  for 
calculating  the  field  in  a  region  where  there  is  a  zonal  structure; 
calculation  according  to  i  mula  (5)  is  simpler  than  a  corresponding 
calculation  by  the  normal-wa  .  e  method  if  the  number  of  normal  waves 
(nondamping  from  r)  is  large. 

Formula  (5)  permits  us  to  calculate  comparatively  easily  the 
acoustic  field  in  the  zones  of  the  refraction  geometric  shadow.  The 
specific  form  of  the  expressions  obtained  as  a  result  of  calculating 
integrals  (6a,  6b)  depends  on  both  the  distribution  of  n2(z)  and 
the  position  of  the  layer  boundaries,  as  well  as  on  the  location 
of  the  sound  source  and  the  point  of  observation. 

Let  us  examine  the  expression  for  the  acoustic  field  in  the 
first  zone  of  the  refraction  geometric  shadow  under  the  condition 
that 

«•(;)•-  c)|  cxp|i>(j(*  .  z)|,  ak„(h  -  *„)* >  ]. 

Let  us  note  that  in  this  case  the  solution  to  equation  (2)  can 
be  expressed  through  special  investigated  functions  so  that 

4-*,7ia|r'>.  i-;.  j r +  4-)]l,’x 

W-i  „  (V"): 

(8) 

** <*s  *) "  -  v  [r  (*  + 1*  +  ) r  (*  ~  i‘  +  4-)r  X 

Xe““  *’  W.Kt(u  lr), 

where 

k'a"’  5-«W  =  2^exp|«(*-i)|. 
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Having  solved  equation  (6a)  by  the  residue  method  we  get 

*«('■.  «.  t,  x„)  exp |-*„(r  -  RJ lm«,|.  ( 9 ) 


where 


'MAoi  *|)  =  °.  0<arg 


An- '  r  2 nr 


4s  iln  | 

(“in-  -J-  j ,ln  ( 

« 

“‘-tJ 

1 

(/  «J  (*«)  —  */  J/  rfl  (#)  ~  *J 

hr  “<*•'>. 

(10a) 


X  exp  +  ikn (r  —  /?oo)  Re  »i  -M  -j-j (I  +  Q,). 

If  (z  -  hQ)  and  (zQ  -  hQ)  are  high  and  l  is  low, 

a  _  l/VT  .  ..■»««■.>  .  t,p  i  T (  or  g°°)lm  I ~  "  <*") I* ) 

lnr  •  (*«.  4  Y »'  ('«>  -  *1  ]/ (*o>  (10b) 

In  equations  (10a)  and  (10b)  zQ  <  h,  z  <  h,  u>0  ^  ■  u(z0,  k,,), 
ui^  =  u(z,  k^),  Iq^I  <<  1,  where 

l 

I #’({)  K-’f’rfS  *“  !  nr.  sln  | /V<'> "f  . 

Ill  h  1  I  fW U )  "  -rj 


-snrcslnse""  *'  j  "l* --y  I  n'M  «f  J ; 

w  (An!  «)  -  -  nrcsln  ^  ■  *'  -Lnlj!!sL_  *'  j  . 

Rm~Rn{2.  »(*„)): 

*)=-  »)  |/I*tt)  •  *T'  ’rfl  I- « 1  |/i5(i)  -.’I-1 2rff. 


(11) 


where 


*lV<S>  I  1  nrcsln'1  *>l  j . 

5,  "  1  r 

f».i  K"(h„)r  |  ■"  |cn;  m|z;  «(A„)|  f*„(r  RJX 

X  llic  x,  •  »<*,)|  -■  (  j’  Rmj  Re  |x,  •  /r(A,i)Pi 

ih  ^n"  I  tlx  2,1 '  v ' 
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In  accordance  with  the  geometric  theory  for  low  values  of 
(z  ~  hg )  the  boundary  between  the  shadow  zone  and  the  first 
illuminated  zone  is  formed  by  a  beam  which  has  touched  the  surface 
z  *  ^0  *  anc*  t*ie  distance  to  the  shadow  zone  boundary  is  equal  to 
R00>  In  this  case  ((3/3K)Rgg)  <  o  and  |Ag^|  decreases  with 
increasing  t.  As  a  result  of  this,  at  r  >’roo  in  formula  (9)  we 
may  restrict  ourselves  to  the  first  terns  of  the  series. 


If  the  value  (z  -  hg)  is  high,  the  shadow  zone  boundary  is 
formed  by  a  caustic.  On  the  caustic  r  *  Rg  ^  z  >  tc  )• 
K0,k  ls  a  root  of  the  equation  (3/3ic)R0(z,  zQ;  k)  -  0.  °  °’k 


5)  >  0,  and  according  to  formula  (10b; 


"00' 

"0,11  lnltially  increases  with  increasing  l  and  only  at 

inf  T  U  4  _ -i  .  .  _  ~ 


In  such  a  case  ( ( 3/3<)Rr 
value  Ag 

sufficiently  high  values  of  l  does  it  begin  to  decrease.  Therefore, 
in  the  vicinity  of  the  shadow  zone  boundary  it  becomes  necessary  to 
take  the  large  number  of  terms  of  the  series  of  equation  (9)  into 
account. 


Now  we  get  an  expression  for  ¥Q(r,  z,  zQ)  convenient  for 
calculation  in  this  case,  and  we  deform  the  integration  contour 
into  a  loop  encompassing  the  poles  of  function  q(K)  in  Integral 
formula  (6a).  On  this  contour 


where 


VE  ■  -y  x 

r.  v  fl’Uo)-**  Y  (!)  —  *» 

x /fc (r.  +  A#(#))t 


^0  (*)  ~r~c  0  (  w  I  A0(»)j  «  1 

at  |  «>  i  >  1 .  |*0r»|»i; 

*.  *0 ;*)  =  */»  +  »„  + ID  is  the  phase. 


(12) 


If  the  shadow  zone  boundary  is  formed  by  the  caustic,  the  phase 
>)>g(r,  z,  Zg)  has  its  saddle  point  kq  in  which  1  -  k2  *  pei$,  where 

*  -  9|  if  r  >  Rgjk;  Kg  =  Kg^k,  if  r  =  Rg  k.  Value  k  Is  the  root 
of  the  equation  ’ 

I  1 1 

V  '><  z'  *<>'  r  -  &>(?.  ?n.  *)  — 0, 
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where  RQ(z,  zQ;  k)  is  set  by  the  equality  of  (11). 


If  z  ■  zQ,  p  ■  exp[a(h  -  z)  -  1; 


*«)  -  v-  ■  -nr.  MnlV'"'  »*  |  r"'*  '~'~lros*  . 

If  0  >  ♦  >  t,  then  Im  Kg  >_  0;  f(Kg)  %  1.  Let  us  plot  the  integration 
contour  in  formula  (12)  through  the  point  k  =  kq  along  the  line  of 
the  most  rapid  decrease  (with  an  increase  of  |k  -  kq|)  of  the  factor 
exp[i<(>0(r,  z,  zQ ;  k)].  Assuming  that  Uq(k0)  -  $0(<e)|  >>  1  and 
I Ag(Kg) |  <<  1,  we  get 

<F„  (r,  z,  z„)  w, , ,  <r.  z,  zj  ■ )-  V  ,,-o  ( (r,  ,Jt  ( 1 3 ) 


where  summation  is  performed  according  to  the  numbers  i  of  the  poles 
touched  in  the  deformation  of  the  contour,  and 


(r,  2,  zB) =-= 


(0(1 

■^I"1  <*„>  ■  4JI,#  i"»  co— ■SJ1'*  x 


X  «P  In  ( r.  r,  i„;  *)  U„  (/,  »)?  x 

J*  -5  *n 


-.  ..f  -'f } 

(I'D 

r  (U  |i'»  ■  ■ 

r  *<>• 

• 

where  v(t)  is  the  Airy  function; 

1  (4°)  [ lb  *)«■--,.]  f~ 

IQo.»l<  1  at  |A, ,(*)„!«  l.  I <•(*„;  «„)|>  I, 

Function  yQ  k(r,  z,  zQ)  at  r  >  R0  k  describes  the  diffraction  from 
the  caustic.  For  an  unbound  medium  lim  't'g(r,  z,  z^)  »  Tq  ^(r,  z*  "g). 

Surface  z  =  hQ  places  limitations  on  the  geometric  dimensions 
of  the  caustic  and  Introduces  perturbations  into  the  field 
diffracted  from  it.  The  series  in  A  in  equation  (13)  considers  this 
perturbation  which  depends  both  on  the  position  and  on  the  reflection 
properties  of  the  upper  boundary.  At  low  values  of  t  we  have 
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*,  z0)=*  Vt't(r,  z,  zJscV^r,  z,  *). 


The  series  In  1  In  formula  (13)  can  also  make  a  noticeable 
contribution  to  VQ(r,  z,  zQ) ,  If  t  >>  1;  In  this  case  fQjlt(r,  z,  zQ)  * 
<  tp  1(r,  z,  Zq )  and  Y0(r,  z,  tQ)  ^  ^Q^Cr,  z»  ZQ)  +  V0,k(r»  z»  20)* 

Let  us  now  examine  the  contribution  to  ¥(r,  z,  zQ)  of  other 
components  which  enter  the  sum  of  expression  (5)  under  the 
condition  hp  ■ 

Let  us  assume  that  n(zQ)  >  n(H)  (the  beam  picture  for  this 
case  Is  depicted  in  Pig.  1).  In  accordance  with  geometric  theory 
the  first  zone  of  the  refraction  shadow  is  limited  on  both  sides  by 
Illuminated  zones  when  n(z)  >  n(H).  If  n(z)  <  n(H),  horizon  z 
intersects  only  the  first  illuminated  zone  and  does  not  pass  through 
the  second  and  subsequent  illuminated  zones,  so  that  the  shadow 
zone  is  unbound  in  the  positive  direction  of  axis  r. 


n/1  jx/a 


Fig.  1.  Beam  picture  (heavy  lines  -  the 
caustic ) . 


Let  us  examine  the  most  interesting  case  when  n(zQ)  >  n(H)  and 
n(z)  >  n(H) ;  we  assume  that  distances  (h  -  zQ)  and  (h  -  z)  are 
sufficiently  small  so  that  the  boundary  between  the  shadow  zone  and 
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the  secondary  illuminated  zone  is  formed  by  the  caustic,  the 
distance  to  which  let  us  designate  by  R1  k.  in  accordance  with 
geometric  theory,  this  caustic  is  the  envelope  of  a  family  of  beams 
whose  equation  has  the  form 


(15) 


where  ■  n(zQ)  cos  a  is  a  real  number  and  a  is  the  angle  of 

beam  emergence  from  the  source. 


•f 

*«••(*.  «•;  «)•=*  !  -  «*r'  ’<*  + « j  ’a 


After  substituting  n(z)  into  the  latter  expression,  we  get 

* f !«'«)  -  ±  amin  . 


The  position  of  the  point  on  the  caustic  is  set  by  the  equation 

where  is  the  root  of  the  equation 

ii *)-°- 


Substituting  equations  (3)  and  (4)  into  formula  (7),  for 
N  "  I i<or<  I  >:>  x>  f(n(H)|  21  1  we  get 


|''£  • 

X  I  I  /‘(K)e  iu"u,,\ || 
k-  F 


(  _  _ 
I  »•'(/»*  *•’  |  n'(t)  ,1  /(,) 

1  'l.l(*>|XPXpi/?,al(r. 

i.i(r,  r,  *>, 


( 16a) 


if  ti/2  <  arg  uH  <  3tt/2  ; 


ft.  1  (r.  ■2rn.1  ■«(;,»)  -  —  ;  «*,-  *k,/l *(| 
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if  f<nrgm< 

/,(*)<-*p[/2»- j.  if  |arr«|<-j- 

m(  )  0(«*  ^  •  *^j;  i-Vi(«)i«  i.  if  |«i> i, 


If  |arg  u„|  <  n/2,  however,  then 


FAr.*,«  111  *7 


V|/*^r  V(*»(l  t  Aj*V) 

I  ?*r  ‘  4  7= - <  / - 

“  K  «»(*»>-«*  y  »»(/)-»» 


=-  X 


X  -  V  p,  (r, 


(16b) 


where 


»  a 

f\.Ar<  *.  4  *„  J  K^fl)  -  t»  d\  4-  kA  Vn'H)  «» A  -  ; 

*  I 

n 

ft.  A'*  *.  *(.:  «)-*/»  4  *#  |  |  *»(|)  ”  «s  rft  4 
«* 

*  _  i 

4-  *.  j  V^f)  »T  <rt  f-  *»  j‘  |  ((T^r <n  - 

*•  *• 

I 

*.  «)-■»  V»4-**J  l/"r(T> :  A  rfl  4 

«  M 

+  *„  j  l/S^T) 

*t  I 

*1  _ _ _ 

*.  «)-*/*+*•  J  K'«,d)-  «,rfi4- 

*  • 

+  2*,  Ij  ^7^*4  *,J  /«»(«)  -«M4-t: 

Iai.|(*)!<  l.where  1«%1>  1,  |*|>1,  |«*|>l,  !*/*!»•• 

Let  us  assume  that  r  <  kJ  then  functions  ^  j(r,  z»  z0*  K^' 

J  ■  1,  2,  3,  H,  5  have  saddle  points  k,  , ,  located  in  the  region  of 

X  f  J 

applicability  of  the  representations  of  formulas  (16a)  and  (16b). 

We  deform  the  integration  contour  of  the  integral  of  equation  (6a) 
described  by  N  ■  lY^r,  z,  zQ)  so  that  it  passes  these  saddle  points 
along  corresponding  saddle-point  contours,  while  outside  the 
essential  saddle-point  vicinities  it  goes  along  the  lines  of  most 
rapid  decrease  P1(r,  z,  zQ;  <).  The  integration  contour  in 
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region  | arg  u^|  ^  it  is  split  into  four  branches  by  the  four 
respective  components  in  expression  (16b);  these  branches  merge  at 
the  origin  of  p(x).  Stemming  from  this,  the  integration  contour 
will  pass  further  through  saddle  point  ie^  Im  ie^  ^  <  0  of  function 
^Cr,  z>  zq»  *)  (Mg.  2).  Saddle  point  < ^  j  satisfies  the  equatior 
(3/3i<)<ti,  .(r,  z,  ic),  which  is  reduced  to  the  form 

1  I J  u 

r  *>.  C7) 


whereupon  function  j(z,  z^;  k )  coincides  with  function  (lb): 

«  ii 

#!./<*.*•;  «)=«]  i«,(«  *rv++,  f  i»’(D 

*»  i. 

+(-lK  '« J|<»‘ («)-»»!  ''M. 


7-2.  3: 


#!.,(*.  *.i  •)“*  j l«*(t) - «»r,,*+ 


(18) 


4- 2.  j  i«>(i)  -  + ( -  i^-'i  i‘  i»»«)  _ 


/« 4,  5. 


Fig.  2.  Plane  <. 

1  —  Point  n(0)  ■  1;  2  —  point  n(H);  3  -  point  n(hg); 

^  -  origin  of  function  ( hQ ;  <);  5  -  pole  of  the 

Integrand  in  formula  (1);  6  —  origin  of  function  p(0; 
7  —  point  k0  ^  and  point  k;  8  -  point  KqJ  9  —  point 

<1  10  -  integration  contour  CQ;  11  -  integration 


FTD-HT-23-439-69 


12 


! 


contour  C^j  12  -  family  of  integration  contours  C,, 

J  *  2,  3,  4,  5.  J 

From  equation  (17)  it  follows  that  0  <  <  n(H)  and  Im  Kj  }  =  0 

for  J  -  2,  3,  5.  Saddle  point  at’r  <  R1>k  lies  in  half-plane 

Im  k  <  0  on  the  line  1  -  l  -  pe4"’,  0  <  *  ^  tt.  Atz-z 

P  *  exp[a(h  -  zQ)]  -  1. 

Correspondingly,  for  R^^Zq,  zq;  k.^)  we  have 

Rt. i (*<>.  *#;  *i, ■)  —  -f  n resin 2e* |  cos_*. 

Calculating  the  integral  of  expression  (6a)  according  to  the  new 
integration  contour  at  r  <  R  we  get 


1f,(r.  *„)=•- V  V,.,(r.  *,  *o), 

/-  • 


J  “  2f  3,  ,  5,  then 


n,(r.«.o--1/  - :x 

V  "? <'«>  «/,/  ( ’  +  Qu,)' 


I'l.  I  (r.  *.  *„)  = 


p  ^  CX|’ ' /Ti' 1  (r'  *•  *-  •>  - 

‘  Tt[ ,J  [33-fli.iU.  *»:  }  • 


I , 1 1  <<  i.  if  lAi,i(<lfl)l  «  i. 
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If  $  <<  1  (small  values  of  R.  .  -  r), 

1  jK 

*u;  *m) -r[^r^i,i(*.  *oi  *).=.1>t](*i. i—  and 
<=s2' 11  *o'J  (/?|,»  —  »,)[^W,1|(*,  *%< 

In  this  case  ^(r,  z,  zQ)  is  reduced  to  the  form 

*.  *o)  = 


-»/» 


[to  *i,  I  <*•  ^  *)■=•!,  ,]  Y (*«)-*’  «  «*  (*)  - 

X*"1, ' 


» 


(19c) 


Formula  (19c)  describes  the  acoustic  field  In  the  region  of  the 
:ic.  It  is  valid  [1]  not  only  with  r  <  R  when  t  >  0,  but 


,k» 

if  t  >  1. 

,  if  t 

i  ('.  z.  z0)  «■--  . 

/ 

/ 

1/  T" 

1 

\  r[  /)» 

X- 

(19a) 


whereupon 


*M<r-  *■  *Di  *1,  |)  -  nrcsln  j  - -  .  iil’  l*  I 

'  t'f 


~  nrcsln  |  1  «JP  l£j*_— *->*">  I  ,  - - - - 

“  I  y'f  e  [4  —  k  Pff'*  -  —  1)2 

+  v  Kpe7*^0 <*-'w~Tj5 . 


It  is  essential  to  note  that  the  use  of  formula  (19c)  at 
t  >>  1  (far  from  the  caustic)  gives  greatly  reduced  values  for  the 
function  ^2>i(r»  z,  zQ)  as  compared  to  the  values  calculated  from 
exact  formula  (19b). 

Let  us  note  that  l'*'ljl(r,  z,  zQ)|  monotonically  diminishes  with 
an  increase  of  (R.  -  r). 

1  ,  K 


1 TD-HT-2  3—^  39—69 


14 


The  integral  of  formula  (6a)  at  N  »  1  also  permits  another 
evaluation,  convenient  at  high  values  of  (R1  k  -  r).  Let  us  set 
the  value  of  r  -  ^  to  which  corresponds  ’ ,  which  satisfies  the 
condition  0  <  $  <<  tt/2.  In  this  case  we  can  plot  the  integration 
contour  in  equation  (6a)  through  ^  without  removing  the  integral 
from  the  convergence  regions  such  that  Im  k  >  Im  on  it.  Then 

we  get  * 

l«Y.(r,  *.  *o)l<  J' I /••,.,(!■.  *.  "f-'Mlx 

r , 

Xj  I  /=■,. i  (r,.  2.  *0;  «)||rf«|, 


whence  it  follows  that  at  r  <  r^ 

I V,. , (r,  *.  ^Xl'tVar,,  20,  *)!e  »•"'  . . I.  (20 

^  In  accordance  with  expression  (5)  Y(r,  z,  zQ )  -  4>0(r,  z,  zQ)  + 
+  ^(r,  z,  zQ).  Comparing  formulas  (6a)  and  (6b)  at  N  »  NQ  ■  1,  we 
see  that  the  integral  expression  for  Y^r,  z,  zQ )  differs  from  the 
corresponding  expression  for  Y-^r,  z,  zQ)  only  by  the  factor 

[1  -  p ( k ) q ( k  )  ]  1  in  the  integrand.  Additional  peculiarities  of  the 
pole  type  Introduced  by  this  factor  into  the  integrand  for 
t'1(r,  z,  zr)  lie  partially  on  the  real  k  axis  at  1  >  k2  >  n2(H), 
and  partially  in  the  regions:  n(H)  £  Re(x)  >  0,  In  it  >  0  and 
0  >  Re  k  >  -n(H),  Im  k  <  0  (see  Fig.  2). 


If  r  <  Rl,k’  the  integration  contour  in  the  integral  of 
formula  (6b)  at  N  ■  1  can  be  deformed  into  exactly  the  same 
integration  path  which  was  used  in  the  calculation  of  Y^(r,  zQ ,  z). 
With  such  transformation  of  the  Integration  contour  the  poles  of 
function  [1  -  p(ic)q(K)  ]-1  are  not  affected  and,  in  addition,  the 
new  contour  runs  in  the  same  part  of  plane  ic  in  which  p(ie)q(K)|  <  1, 
if  | V(k) |  <  1.  After  deforming  the  integration  contour  into 
contour  Cj  by  the  method  shown  and  after  estimating  the  value  of  the 
Integral  we  get 


(21) 
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where 


In  other  words , 

Q,. ,  l<  !  exp [•*.*  >  r)1'* (£  /?,. ,  («. *;  J-«j -1  j  •;  (  22a) 


I  Qi.  / !  c>  7~.  •  JL  2*  3-  5. 


Estimation  of  equations  (22b)  Is  convenient  if  | V( >e .  )|  <<  l. 
expression  for  ^(r,  z,  zQ) ,  J  -  2,  3,  Hp  5  can  also  be  obta 

*•  *«>=-  22  x.  x„). 


In  the  form 


(22b) 

The 
obtained 


(23) 


where 


Tx./l'.  *•  *<■) 


'I /  - ;•  X 

r  f(rf.  wa. 

y  *»■  »*.<>  •  j  (AT  I)  '{.j 


(?n> 


Here 


T*.,(r.  t.  X,;  *)*=>ti,#(r,  *,  «)  f  2(Af  «>; 

H 

*"•/<*•  *»;  «)-=#,.,<*,  «)  f-2(W  l)«  j  |n*(») 

CN,J  ls  the  root  of  the  equation  r  ■  R„  ,(z,  z_;  tc). 
*  •’•JO 


(25) 


Formula  (23)  is  obtained  by  replacing  [1  -  pU)q(ic) J"1  by 
J_i[p(K)q(<)]  in  the  integral  describing  ^  j(r,  z,  zQ)  and  the 

tiext  obtained  integration  of  the  series  obtained.  The  integration 
proceeds  along  the  contour  obtained  by  deforming  the  original 
integration  path  into  the  line  of  most  rapid  decrease  exp[i<$  (p  ■/. 
<)]  with  an  increase  of  \k  -  k  | .  This  line  passes  through 
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saddle  point  *N  j  and  ends  at  the  origin  of  function  p(ie).  Function 

1>  ,  (r,  z,  zn)  describes  a  wave  corresponding  to  a  geometric  beam 

N  f  j  U 

reflected  N  times  from  the  boundary  z  ■  H. 

The  expression  in  the  right  side  of  (22a)  is  considerably  less 
than  unity,  even^with  small  values  of  ^  -  r).  It  is  easy  to  be 
convinced  that  | Qx  <<  1  is  also  much  less  than  1  at  r  ■ 

Having  calculated  the  values  of  V^Cr,  z,  Zq)  and  4^(r,  z,  Zq) 
by  the  method  shown  above  we  can  now  find  the  complete  field 
Y(r,  z,  zQ),  from  formula  (5),  in  the  geometric  refraction  shadow 
zone.  In  Fig.  3  there  is  shown  the  dependence  of  function 
»(r,  zQ)  on  r,  calculated  for  the  particular  case:  a  *  u-10  1/m, 

n  -  z  -  100  m,  H  -  z  ■  800  m  under  the  assumption  that  pg  ■  p;  under 
this  condition  the  contribution  of  the  waves  reflected  from  the 
horizon  z  ■  H  to  field  ¥(r,  z,  z Q)  is  negligibly  small.  In  the 
graph  of  Fig.  3,  the  points  mark  the  values  of  Y(r,  z,  zQ), 
calculated  by  the  normal-wave  method.  The  results  obtained  by  both 
methods  give  satisfactory  agreement. 


W««= 


Fig.  3.  Graph  of  the  dependence 
of  function  V(r,  z,  zQ)  on  r. 

1  -  Values  of  |f|  calculated  by 
the  normal-wave  method;  2  -  shadow- 
zone  boundaries. 

Let  us  analyze  the  case  when  n(z)  <  n(H),  u(H,  n(z))  >  1.  Here 

'V 

the  integral  of  formula  (6b),  which  describes  ^(r,  z,  Zq),  may 
conveniently  be  estimated  by  the  stationary-phase  method. 

Let  us  formally  expand  the  integrand  in  formula  (6b)  into  a 
series  of  powers  of  (p(<)q(<)).  The  terms  of  this  series  have  real 
stationary-phase  points  only  at  k  t.  [0,  n(z)],  for  which 
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.  w  R 

t  \  {r,  t,  t„:  *.)  \i  y*  _ 

r  -/> (»)»'(*)  ;  *,  *0;  *)l/»(*)^(«)iAr“l, 


(?6) 


while  z,  Zq  ;  k  )  satisfy  equation  (l6b).  The  essential 

stationary-phase  points  which  correspond  to  a  term  of  the 

series  in  the  right  side  of  formula  (26)  are  roots  of  equation  (25). 

Let  us  plot  the  integration  contour  so  that  Im  k  =  lim  (+E)  on 
it.  Then  e-° 


(r'  *«  *">  ■  ^  *•  *o)+t  (r ,  «,  *0). 


(27) 


Here  ^N, J ( r*  z*  V  corresponds  to  the  contribution  of  stationary- 
phase  points  of  a  term  of  the  series  of  formula  (26). 

If  equation  (25)  has  a  unique  root  <N  then  ¥  (r,  z,  8  )  is 

described  by  equality  (24).  0 

But  RN,j(z*  z0i  K)*  where  J  =  3,  5,  as  a  function  of  real  values 
of  has  a  maximum  at  point  k  -  <NJ>kJ  which  corresponds  to  the 
caustic  of  the  beams  reflected  from’the  horizon  z  «  H.  This  caustic 
lies  completely  in  the  region  of  values  of  z  for  which  n(z)  <  n(H). 
Magnitude  <NJ(]<  is  the  root  of  equation  (z,  Zq;  k)  =  0.  The 

distance  from  the  sound  source  to  the  caustic  along  the  horizontal 
is  equal  to 

Rn,  i,  t  —  7?/v,  j(z,  z0;  **,  j > 


If  r  is  close  to  R, 


and  r  >  RN,j,k‘  then 


. ! ;  aS***W.  »  *  V<v./)e»P|/fA-.  ;(«••*.  ».A.  y  ,) | 

<V."  #,v. J  |/  "JUu) 


(26a) 


where 
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M  , 

Function  l  ^(r,  z,  z„)  is  the  sum  of  normal  waves  corresponding  to 
m=  1 

real  poles  [n(h),  n ( H ) ] .  A  normal  wave  is  defined  as  the  residue 
in  pole  Km  of  the  Integrand  in  formula  (1). 

With  the  given  form  or  n2(z),  ?m(r,  z,  zQ )  is  described  by  the 
formulas  given  in  worker].  The  function  X?m(r,  z,  z  )  gives  the 
basic  contribution  to  ^(r,  z,  zQ),  if  |v|  <<  1  or  |V|  <  1,  r  ■*  », 
and  characterizes  the  diffraction  "exposure"  of  the  geometric 
refraction  shadow  zone  stretching  limitlessly  in  the  positive 
direction  of  the  r  axis  on  horizon  z  under  consideration. 

Although  at  n(z)  <  n(H)  horizon  z  does  not  pass  through  the 
secondary  and  subsequent  illuminated  zones,  the  resultant  field 
fCr,  z,  2q )  will  nevertheless  be  characterized  on  it  by  alternating 
maxima  and  minima.  The  picture  is  externally  similar  to  that  which 
take.'  place  at  n(z)  >  n(H).  But  in  the  case  in  question  the 
appearance  of  maximum  |f(r,  z,  zQ ) j  is  caused  by  the  approximation 
to  the  caustic  of  the  beams  reflected  from  the  ground,  instead  of 
the  input  to  the  illuminated  zone.  Such  a  false  "illuminated  zone" 
differs  from  the  true  one  in  that  maximum  ^(r,  z,  zQ)|  in  this 
zone  depends  on  the  frequency  of  the  sound,  even  with  the  consideration 
of  absorption  in  the  medium. 

The  acoustic  field  in  the  channel  is  represented  in  the  form 
of  a  finite  number  of  components,  each  of  which  is  described  by  a 
converging  integral.  This  representation  is  especially  convenient 
in  the  region  where  the  field  has  a  zonal  structure,  since  it 
permits  us  to  calculate  the  field  here  faster  and  simpler  than  by 
the  normal-wave  method  ordinarily  used  for  this  purpose. 

The  field  was  calculated  in  the  refraction  geometric  shadow 
zone  in  which  n“(z)  =  2  exp[a(h  -  z)]  -  exp[2a(h  -  z)].  It  Is 
important  to  note  that  the  obtained  formulas  are  also  valid  with 
another  similar  dependence  of  n2(z),  since  the  specific  form  of  n2(z), 
as  a  rule,  is  not  taken  into  account  in  the  derivation  of  the 
formulas . 
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